LTH 795, DESY 08-092 (July 2008) 



00 
O 

o 



oo 

Oh 1 



> 
00 

m 

(N 

o 

oo 
O 



X 



Towards the NNLO evolution of polarised parton distributions* 
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We report on the first calculation of the structure function g 1 in polarised deep-inelastic scattering to the 
third order in massless perturbative QCD. The calculation follows the dispersive approach already used for the 
corresponding unpolarised cases of JVl, but additionally involves higher tensor integrals and the Dirac matrix 
7 5 in fl / 4 dimensions. Our results confirm all known two-loop expressions including the coefficient functions 
of Zijlstra and van Neerven not independently verified before. At three loops we extract the helicity-difference 
next-to-next-to- leading order (NNLO) quark-quark and gluon-quark splitting functions AP qq and AP qg . The 
results exhibit interesting features concerning sum rules and the momentum-fraction limits x — > 1 and x — > 0. 



1. Introduction 

Two decades after the seminal EMC measure- 
ment pQ set off the so-called spin crisis, the inter- 
nal spin structure of the nucleon remains a very 
active field of research. By now the measurements 
of polarised (spin-dependent) deep-inelastic scat- 
tering (DIS) - the very process first measured by 
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gP(x,Q 2 ) + C(x) 



= 0.0063 (+7.5) 

k = 0.0141 (+6.2) 
x = 0.0245 (+5.2) 
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0.0346 (+4.5) 
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x = 0.0775 (+3.5) 
0.122 ( + 3.0) 
= 0.173 (+2.5) 
= 0.245 ( + 2.0) 
x = 0.346 (+1.5) 
x = 0.490 (+1.0) 
x = 0.735 (+0.5) 
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Figure 1. The end-2006 world data on the proton 
structure function g 1 with some NLO QCD fits [5] 
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EMC - have reached a high accuracy but are still 
restricted to relatively small scales Q 2 , see Fig. 1. 

Dramatic further improvements, in particular 
of the kinematic coverage, can be expected from 
the planned Electron- Ion Collider, sec Rcf. [3]. 
But, as rather low values of Q 2 imply rather large 
QCD corrections, even analyses of available data 
will profit from calculations beyond the present 
next-to-leading order (NLO) approximation. 

2. Polarised DIS in perturbative QCD 

The kinematics and perturbative-QCD factor- 
ization of inclusive DIS are recalled in Fig. 2. 

Scale, Bjorken variable 

Q = -Q 
x = Q 2 /(2p-q) 

n (p) (fT \ Lowest order : x = £ 

Figure 2. DIS in the QCD-improved parton model 
The measurements of the difference a^^ — a^^- 

e n en 

of the hclicity-dcpcndcnt cross sections yield the 
polarised structure function of the nucleon n. 
In perturbative QCD, neglecting contributions 
suppressed by 1/Q 2 , this quantity is given by 




5i"(x,Q 2 ) 



A/£(£,M 2 ) 
(1) 



1 
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where, as in Eq. |J2J) below, the sum over the 
parton species k is understood. c 9r k represents 
the corresponding coefficient functions (mass- 
factorized partonic cross sections) for g 1 , a s de- 
notes the strong coupling, and \i = 0(Q) stands 
for the renormalization and mass-factorization 
scale. The helicity-difference parton distributions 
Afi = fi—> — obey the evolution equations 

A/ 4 ( M 2 ) = [AP lk (a s (^)) ® A/,( M 2 )] (2) 

where ® denotes the Mellin convolution written 
out in Eq. (fT]). The initial conditions for Eq. @ 
are, of course, incalculable in perturbative QCD. 
Present lattice techniques can provide, at best, 
only a few low Mellin moments. Hence predic- 
tions for, e.g., processes at RHIC, can be obtained 
only via fit-analyses of reference obscrvablcs such 
as g 1 and the universality of the parton densities. 
For the present state-of-the-art see, e.g., Rcf. [3]. 

The perturbative quantities entering a general 
hard-scattering observable a are thus the splitting 
functions P and the coefficient functions c a , 

P = a 5 P(°) + a 2 P« + a 3 s P^+ ... 

c a = a?-\<g>+a B £>+ci<& + ...]. (3) 



Being the first serious approximation, the NLO 
cannot provide a sound estimate of the error from 
the truncation of the perturbation series. The 
NNLO requires the two-loop coefficient functions 
for all observables under consideration and the 
universal three-loop splitting functions. The lat- 

(2) 

ter quantities AP> k for the polarised evolution 
([2]) are the main subject of the present research. 

3. Previous second-order calculations 

Before we turn to our new third-order calcu- 
lations, it is useful to briefly recall the previous 
two-loop results for the coefficient functions for 
g x and the splitting functions AP ik . All those 
calculations have been carried out in dimensional 
regularization. Hence the helicity-difference pro- 
jection of initial quarks introduces the issue of 7 5 
in D 7^ 4 dimensions into the problem. 



The perturbative corrections for the polarised 
structure function g 1 to order a 2 have been calcu- 
lated 15 years ago by Zijlstra and van Neerven [5]. 
This x-space calculation, carried out using the so- 
called Larin scheme [6] for j 5 (see below) , yielded 

(2) 

the NNLO coefficient functions c , and well as 

9i,q/g 

two of the four polarised NLO splitting functions: 
AP^ and Apjg } . 

The complete set of NLO splitting functions 
APjp was first computed two years later by Mer- 
tig and van Neerven using the operator-product 
expansion [7J. Here 7 5 was treated using the 
reading-point method of Ref. [5] ■ The results were 
independently verified already a few months later 
by Vogelsang in an axial-gauge calculation. 
This calculation was performed using the origi- 
nal HVBM prescription for 7 5 [TU] but included 
checks with the other schemes mentioned above. 

All three theoretically (but not computation- 
ally) identical prescriptions for 7 5 require an ad- 
ditional renormalization or factorization of the 
loop results in Refs. [5|7|9j . The complete set 
of the corresponding expressions for g 1 to NNLO 
has been derived in Ref. [TT] ten years ago. We 
will return to this issue below. 

4. The three-loop computation via forward 
Compton amplitudes 

The present calculation follows the approach 
used before in our computation of the unpolariscd 
splitting functions [12] (see also Rcf. [13]) and 
the coefficient functions for F^.l [33] to order o? s . 
Thus we employ the optical theorem 




to transform the 7*/ total cross sections into for- 
ward Compton amplitudes. Due to a dispersion 
relation in the partonic Bjorken variable x, the 
coefficient of (2p-q) N provides the N-ih moment 

F a N = jld XX N -'F a {x) 

of the partonic structure function F a selected by 
a suitable projection of the 7*-amputatcd graphs. 
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In particular, the projection of the partonic 
tensor W^" on g x in D dimensions reads [5] 



& = 2[(D-2)(Z>-3)(p.g)]- 



(4) 



using the Schoonschip notation for contracted 
indices. The N !_1 LO splitting functions AP^ 

1M 



and the N'LO coefficient functions c„ \ can be 
extracted from the l/e poles (D = 4 — 2;) and 
e° parts, respectively, of the corresponding Z-loop 
results for 7*/ -» 7*/, cf. Refs |l2ll3ll4j . 

Treatment of the integrals 

In order to illustrate the computation of inte- 
grals required to evaluate the forward Compton 
amplitudes, let us consider one of the hundreds of 
7*<7 three-loop diagrams shown here together with 
a pictorial representation of its momentum flow: 



For the latter one disregards the external parton 
lines and draws the remaining self-energy type 
diagram, the topology of which is denoted follow- 
ing the notation of Refs. [THj . Our example is a 
ladder (LA) diagram. The (partly additional) de- 
nominators carrying the parton momentum p are 
then indicated by the fat (in the coloured version: 
red) lines. Here p runs, after turning the diagram 
upside-down, through the lines 1, 2 and 3. Thus 
the example is assigned the subtopology LA13. 

According to our above discussion, we need an- 
alytic expressions for the (dimcnsionless) coeffi- 
cients I{N) of (2p -q) N /Q 2a . One might try to 
obtain I(N) by brute force, Taylor-expanding the 
denominators with p and working out the sums. 
It turned out that such a strategy, in general, 
does not work. Instead, we employed identities 
based on integration by parts, scaling arguments 
and form-factor decompositions (see Sect. 2 of 
Ref. [12] 1) to successively simplify the integrals. 

The LA13 integrals, e.g., can be simplified by 
applying p^d/dq^ both inside and outside the 
integral. For the scalar integral with unit denom- 
inators this yields 



<n> 



iV + 3 + 3e 2p-q 
N+2 q 2 



<n> 



N+2 



<n> 



(5) 



Here the LA13 integral occurs twice, once with a 
prefactor 2p -q. Hence Eq. (5) represents a differ- 
ence equation (here of order n = 1) expressing its 
coefficient I(N) in terms of that of a LA12 inte- 
gral with an enhanced denominator in the 3-linc, 

a (N)I(N) - ... - a n (N)I(N-n) = G(N) . (6) 

First-order recursion relations like Eq. (5) can be 
reduced to a sum. Higher-order recursions (we 
needed equations up to n = 4) can be solved by 
inserting a suitable ansatz into Eq. ([5]). Both pro- 
cedures exploit the fact that all integrals required 
for the computation of the splitting functions can 
be expressed in terms of harmonic sums |16| . 

Despite being uncharacteristically simple in 
both derivation and size, Eq. (5) illustrates the 
strict hierarchy of subtopologies in our procedure. 
Our LA13 example can be evaluated only once 
the LA12 integral in Eq. (5) is known. This in- 
tegral requires the so-called basic building blocks 
(with only one p-dependent denominator) LAn 
and LA22 together with other integrals of simpler 
topologies where a non-p denominator has been 
removed. Also those integrals need to be evalu- 
ated in terms of yet simpler cases, and so on. 

During the first half of this decade, reduction 
chains for all subtopologies were constructed and 
coded in Form [17]. This and the computation 
of all integrals (and their much more numerous 
sub-integrals) entering the diagrams for the unpo- 
larized cases [12113114] took years of both human 
and computing resources. It would not have been 
possible to get through without extensive tabula- 
tion of intermediate results. By 2005, a database 
had been accumulated of more than 100 000 inte- 
grals requiring about 3.5 GBytes of disk space. 

However, as we will see below, this database 
is still not large enough to cover all integrals re- 
quired for the three-loop calculation of the spin- 
dependent structure function g 1 . 
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Numerators for a non-planar diagram 

One of the more difficult three-loop diagrams 
for the 7*9 Compton amplitudes is shown below 
in fully amputated form: 



2 




5 



l{il,P7, i 2,P4,! J ->P3> i 3>P8> i 2,P6' l 'iPli i l>P+P2> i 3) 

Pi ■ ■ -pi {P+P 2 ) 2 (.P + P7) 2 

Here the numerator has been written in the no- 
tation used by Form [TTj, and the numbers in 
the graph correspond to the line numbering con- 
vention in Mincer [TS]. Thus this diagram, with 
the quark entering in line 2 and leaving in line 7, 
is of the subtopology NO27. 

In order to determine the unpolariscd splitting 
function, it is sufficient to evaluate (7) ■ j{p) g„ v . 
The resulting trace leads to numerators 

(P2-P) k2 (p 3 -P) k3 (pr<l) k9 with fc 2 + fc 3 + fc 9 <3 

and the denominator l/f>| cancels in all terms, 
which already reduces the number of required re- 
duction steps by one. 

The disentangling of the contributions to the 
coefficient functions for P 2 and Fl required the 
additional calculation of (7) • j(p) p^Pv / '(p ■ q) 2 
which includes higher numerators, 

O2 -p) k2 (P3 -p) k3 (P2-i) k9 with k 2 + k 3 + k 9 < 4 . 

Again no denominator 1 jp\ remains in any term. 

The determination of the spin splitting func- 
tion and the coefficient function for g 1 requires 
the evaluation of (7) -"/(p, 5) e^ pq /(p-q) with the 
7-term arising from the quark hclicity-diffcrcncc 
and the £-tensor from the g 1 projection in Eq. (J3J). 
Now one encounters yet higher numerators 

(P2 -P) k2 (P3 -P) k3 (P2-l) k9 With fc 2 + k 3 + kg < 5 

and the denominators 1 jp\ remain in terms up to 
the highest numerators. Thus, compared to the 
previous calculations of F2 and Fl [12113114] , two 
additional reduction steps (often leading to addi- 
tional denominator enhancements, recall Eq. (5)) 
are required in the polarised case. 



The situation is similar in other difficult sub- 
topologies including NO12 , BEgs and LA17 (some 
of which, unfortunately, are far from being fully 
automated), leading to a large number of new 
(sub-) integrals beyond the 2005 database. 

Using codes and experience from the previous 
unpolarizcd calculations, we were able to compute 
the required integrals within a couple of months. 
This would have been impossible without the 
possibility of fixed- checks at all intermediate 
stages facilitated by the Mincer program |15j . 

Treatment of 7 5 and two-loop checks 

We employ the Larin prescription for the quark 
helicity-difference projector, i.e., we use [6] 

1 

i>l5,L = ^^vpltilvlp ■ (8) 

Hence, after coupling-constant renormalization 
and mass-factorization, we need to perform a 
scheme transformation Z which nullifies the first 
moment of the non-singlet (ns) quark-quark split- 
ting function in accordance with the conservation 
of the non-singlet axial vector current [719] , 

9i = Cg i>L Af L 

= c ai>L Z- 1 ZAf L = c g - s Af- s . (9) 

The general form of the corresponding transfor- 
mation of the splitting functions to NNLO in the 
flavour-singlet case, where Eq. © has to be read 
as a matrix equation, can be found, e.g., in Sect. 2 
of Ref . [18] . For the special case 

Z - - 5 - +a z^S + a 2 (z {2} + z {2) )S (10) 

— <->lj 1 u s^qq Uqq i "s V^(jg,ns ' ^qq,psj u qq ) 

the (x-space) transformations of the NLO and 
NNLO parton-quark splitting functions read 

8[AP$] = -faff , 

5[AP«] = 4>®AP(°) (11) 
and 

6[AP$] = M ( Zqq ^ 2 - 2z$] - (3 lZqq \ 

6[AP$] = 4 2 >ApW + zW0AP« L .(12) 

The non-singlet entries in Eq. (fT0|) can be in- 
ferred from the fact that the MS coefficient func- 
tions for g 1 and the structure function F^ +v in 
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neutrino-nucleon DIS arc identical to two loops. 
The functions thus determined agree with the re- 
sults of the direct calculation in Ref. [11] . At the 
moment the latter article is the only source for the 

(2) 

NNLO pure-singlet transformation term Zqq^ s . 

After transforming to the cc-space expressions 
in terms of harmonic polylogarithms |19] and per- 
forming the transformation ([9]) - (jTTJ) . our two- 
loop expressions fully agree with the previous re- 
sults for AP„l and discussed in Section 3. 

q 91 ' (2) 

Incidentally, we note that the change of c gitPS in 
the second erratum of [5] is exactly the scheme 
transformation induced by z qq ] ps of Ref. [llj . 

5. Checks and features of the 3-loop results 

A strong check of the third-order computations 
is provided by the e~ 3 and e~ 2 pole terms of the 
unfactorized partonic structure functions. These 
terms have to be of a specific form in terms of 
the lower-order splitting functions and coefficient 
functions, including the NLO quantities APgf of 
Rcfs. |7I9| inaccessible to a two-loop calculation of 
photon-exchange DIS. Our results pass this test. 

We find that AP„ ( s } in Eq. @ is equal to the 

(2) — 

unpolarizcd non-singlet splitting function Pn S 
after the scheme transformation (fl~2]) adjusting 

Cgjns to the (y + v) coefficient function c F [20] , 

^ (2) 

This guarantees that the first moment of AP ns 
vanishes as required by the conservation of the 
non-singlet axial vector current. 

Also the first moment of the NNLO gluon- 
quark splitting function vanishes 

AP^\N=l) = (13) 

in agreement with the 1990 result of Altarelli and 
Lampe [21]. For the pure singlet the relation 

AP p ( s n) (A^=l) = -2n f AP^- 1] (N = l) (14) 

also holds at NNLO (n = 2). Eqs. ([13]) and ([14]) 
suggest the possibility of a higher-order generali- 
sation which would bring us closer to a complete 
NNNLO description of g 1 at N = 1. 

Turning to large x, we note that the helicity flip 
is suppressed by two powers of (1— x) as x —>■ 1, 

P${x)-ApW{x) ~ (l-xf-P^(x) (15) 



(the corresponding relation for AP q q ' is obvious 
as only the suppressed pure-singlet parts differ), 
where P qg denotes the unpolarized (flip + non- 
flip) splitting function. This is in line with the 
pattern advocated by Brodsky et al. [22] for the 
perturbative and non-perturbative contributions 
to the helicity-dependent parton distributions. 
It is interesting to note that the same suppression 
holds for 7 of the 8 LO and NLO splitting func- 
tions AP^' X \ the curious exception being APgq^ 
(we are not aware of any earlier discussion of this). 
A simple scheme transformation (zgq oc APgq ) 
exists which would 'repair' this exception. 

Since it is identical to the unpolarized quantity 

(2)— (2) 

P ns , the leading small- a; logarithm of AP ns 
agrees with the prediction of the small- a; resum- 
mation [23]. The situation is less straightforward 
in the singlet sector where we find 

AP £\* X = -5Cjn f -IC F n f (C A -n f ) 

(Ca and Cp are the usual group factors of QCD). 
The first of Eqs. ([16]) agrees with the 1996 predic- 
tion of Bliimlein and Vogt [24] based on the sin- 
glet resummation by Bartels et al. [25]. However, 
the corresponding coefficients of CaCftt-j: and 
Cp-rij: of AP q g 2 ^ are not the same. Agreement is 
restored in the supersymmctric limit. 

We have to leave the clarification of this point 
to further studies. However, we should note here 
that, unlike in the non-singlet case, the leading 
small-x logs of the singlet splitting functions are 
not scheme invariant. This can be seen from 
the second lines of Eqs. (fTT]) and (fT2"]) : for z^ q oc 
Cg™q, for example, <$[AP q ^] will include leading- 
log contributions (recall that ln fc x ~ 1/A^ fe+1 in 
A^-space). Hence one cannot conclude now that 
the above difference constitutes an inconsistency. 

A final remark is in order on the 7 5 -induccd 

transformation to the MS factorization scheme. 
(2) 

An additional term z gq can be expected from the 
use of, e.g., the Larin scheme in a calculation of an 

(2) 

observable suitable to access AP g q (see below). 
This presently unknown function would add an 
additional term —Zgq <S> AP q g' to the first line of 
Eq. (fT2"]) without affecting the result for AP q ^. 
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Fig. 3 finally illustrates the pcrturbative expan- 
sion of the hclicity-diffcrcncc gluon-quark split- 
ting function to the ncxt-to-next-to-lcading order. 
Here we have employed the code of Ref. [25] for 
the numerical evaluation of the harmonic poly- 
logarithms [19] up to weight four. 




-3 -2 -1 

10 10 10 1 

X 



Figure 3. The LO, NLO, and NNLO approxima- 
tions for the splitting function AP qg at values of 
a s and the number of light flavours relevant to 
present data on the spin structure function g 1 . 

The third-order corrections are 15% or less over 
the wide region 0.005 < x < 0.9 even for this rel- 
atively large value of a s . The figure also shows 
that, as for all other three- loop quantities we 
have calculated |12|13|14j (see also Refs. [27] for 
the time-like case), the leading small- a; logarithm 
docs not provide a useful approximation to the 
complete function at any relevant values of x. 



6. Summary and Outlook 

We have calculated the three-loop QCD cor- 
rections to spin-dependent electromagnetic deep- 
inelastic scattering and extracted the upper row 

(2) 

AP f (x), f = q,g, of the matrix of the third- 
order (NNLO) spin splitting functions. 

Due to the presence of difficult tensor integrals 
beyond those encountered before [12113114] . the 
computation of the corresponding diagrams pro- 
vided yet another example of the 'no free lunch' 
theorem. The D ^ 4 continuation of 7 5 has been 
dealt with using the Larin prescription [6] and, 
as in previous two- loop calculations [719] . a fi- 
nal scheme transformation to MS factorization. 

(2) 

The result for AP qg is final also in this respect, 
while the pure-singlet quark splitting function 

(2) 

AP ps may require an additional transformation. 

Our calculation also includes the coefficient 
functions for the spin structure function g 1 up 
to the NNNLO (again in MS up to a scheme 
transformation in the pure-singlet sector). At two 
loops we confirm the previous results for the 

(2) 

NNLO coefficient functions c g f of Ref. [S] which 
were not independently checked before. Return- 
ing to the NNLO splitting functions, we note that 
our results exhibit interesting features at the mo- 
ment N = 1 and the momentum-fraction limits 
x — > 1 and x — > which merit further studies. 

(2) 

We have started work on the lower row AP g l f , 
/ = OS 9 °f the splitting-function matrix. Our 
method is a generalisation of the scalar-exchange 
approach used in the unpolarized case [12113] . i.e., 
we now consider DIS including also the exchange 
of a pseudoscalar (Higgs in the heavy-top limit) % 
coupling to G^ U G£ U . This adds an operator mix- 
ing to the problem [28] . Moreover, there are many 
more diagrams in this sector, e.g., a preliminary 
(non-optimised) database has about 30000 g<figx 
diagrams where 4> is the scalar Higgs used before. 

Despite the expected impact of improved tools 
such a the new multi-threaded version of Form 
[29] . it will thus take some time before all NNLO 
spin splitting functions are known. In the mean- 

(2) 

time, the full TV- and x-space results for AP q j ; 
(including Form and Fortran routines) will be 
presented in a forthcoming journal publication. 
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